Effective a Model Formulation for Two Interacting Electrons in a Disordered Metal 
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We derive an analytical theory for two interacting electrons in a d-dimensional random potential. 
Our treatment is based on an effective random matrix Hamiltonian. After mapping the problem 
on a nonlinear a model, we exploit similarities with the theory of disordered metals to identify a 
scaling parameter, investigate the level correlation function, and study the transport properties of 
the system. In agreement with recent numerical work we find that pair propagation is subdiffusive 
and that the pair size grows logarithmically with time. 
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In the theory of Anderson localization one is accus- 
tomed to characterize localized single particle states by 
their typical extension. This so-called localization length 
is also believed to be the relevant length scale separat- 
ing metallic and insulating behavior. However, it has 
become clear recently that two interacting electrons can 
form a pair that propagates up to a scale L2, far beyond 
the one-particle localization length L\. This effect, orig- 
inally proposed by Shepeleyansky jjj and subsequently 
investigated, made more precise and generalized in a se- 
ries of papers has added an important aspect to 
the theory of localization. Inhowfar this phenomenon is 
amenable to experimental verification is presently under 
active consideration in the community. The theoretical 
information accumulated so far is, however, mainly of ei- 
ther qualitative or numerical nature. It is the purpose 
of this letter to develop an analytical theory of the effect 
for length scales L > L-y, where the interaction-assisted 
pairs are well defined. 

Starting from an effective Hamiltonian we derive a non- 
linear a model describing the problem in arbitrary di- 
mensions. This enables us to address a number of inter- 
esting issues. We identify a certain scaling parameter, 
namely an effective "pair conductivity" <J e ffi which co- 
incides in one dimension with the pair localization length 
L2 ~ & e ff ~ L\. In the perturbative regime, the a model 
gives rise to subdiffusive dynamics, where both the diffu- 
sion constant T> e ff(uj) of the pairs and the local density 
of diffusing pair states v e ff{u>) depend on frequency. As 
a consequence, we find that (for a certain time scale) the 
pair size grows logarithmically with time and that the dif- 
fusion is supressed by a similar logarithmic factor. These 
results have some very interesting relations with recent 
numerical findings for two "interacting" kicked ro- 
tators. Finally, we identify two regimes with different 
spectral statistics separated by a pair Thouless energy 

The first step is to motivate the effective random ma- 
trix Hamiltonian (ERMH) that serves as our starting 
point. We consider two interacting electrons on a lattice 
in d dimensions with site diagonal disorder and nearest 



neighbor hopping elements. We assume that all one- 
electron eigenstates are localized with a typical localiza- 
tion length Li > 1 (all length scales are measured in 
units of the lattice spacing). For d = 1,2 this condi- 
tion implies sufficiently weak disorder while for d > 2 
the disorder strength should be slightly above the crit- 
ical disorder W c . The localized one electron states are 
denoted by ip p (r), where the index p refers to the "cen- 
ter" of the state. The Hamiltonian, expressed in the basis 
of symmetrized (i.e. both electrons have opposite spin) 
one-electron product states \p\p2) with {r\T2\p\P2) — 
[<P P i(ri)<P P2 (r2) + ip Pl (r 2 )fp 2 (ri)]/V2, reads 
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Here, the Q 



P1P2, P3P4 



are the interaction matrix elements 
expressed in the basis of the product states \p1p2 > and 
e p is the one-electron energy of the localized state <p P (r). 
We emphasize that we consider the dynamics of the pairs 
in the basis of the above product states. Only for scales 
L > L\ can this be easily interpreted as motion in real 
space. For a local Hubbard interaction with interaction 
strength U we have tjj: 



(2) 



These interaction matrix elements become exponentially 
small whenever any two of the "positions" pj differ by 
much more than L\. Shepelyansky mapped [Q (for 
d = 1) the Hamiltonian (jl]) on a random band ma- 
trix with a superimposed, strongly fluctuating diagonal 
matrix. To do so he employed two crucial approxima- 
tions: First he neglected all badly coupled pair states 
(those with \pi — p 2 | > -^l) an d second he assumed 
the remaining coupling elements (^|) to be independent, 
equally distributed random variables. With the ansatz 
</?p( r ) ~ e - ' r ~ p !/ Ll a p (r) / L^ 2 , where a p (r) is a random 
variable with (a p (r)} = and (a p (r)a,p(f)) = S pp S r f, 
and using the central limit theorem one gets the esti- 
mate §| (Q 2 piP2tP3Pi ) - U 2 L^ 3d for the well-coupled 
states. Our ERMH defined below only relies on the sec- 
ond of the above assumptions: All pair states are taken 
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into account, but the coupling matrix elements (g) will 
still be taken to be independent Gaussian variables. This 
latter point is indeed quite a serious simplification since 
the correlations between different coupling matrix ele- 
ments are neglected. We believe, however, that essential 
physics can be learnt even without paying attention to 
this refinement. In particular, the strongly fluctuating 
diagonal elements tend to eliminate the effect of those 
correlations. 

Let us consider as configuration space a 2d- 
dimensional lattice with two coordinate vectors R = 
Pi + 92 and j = (pi — p2)/2 corresponding to twice the 
center of mass and half the distance of the two elec- 
trons, respectively. Our ERMH, Ji = f) + consists 
of a strongly fluctuating diagonal part i) with entries r, R 
and of an interaction induced coupling matrix The 
?7g correspond to the one-electron energies e Pl + e P2 in 
(M. We take them to be independent random variables 
with the distribution function po(v)- We typically have 
p (rj) ~ l/(2W b ) for |r?| < W b , where W b is the band- 
width of the disordered one-electron Hamiltonian. The 
matrix elements of £ are independent Gaussian random 
variables with zero mean and variance 



( c fj - 

\*rr 



.a(\R-R\)v(j)v(j) 



Here a(\R — R\) and v(j) are smooth functions decay- 
ing exponentially on the scale L\. We need not specify 
their particular form, it is sufficient to know their typical 
behavior 

' U 2 L^ d ,|#|<£i, 



1 ,\j\<L 1: 



(4) 



which is justified by (||) and the exponential decrease 
(on the scale Li) of the localized eigenfunctions. The 
function a(\R — R\) describes how the coupling strength 
decreases exponentially if the distance of the centers of 
masse increases, wheras v(j) describes how the increasing 
size of the pair states reduces the coupling. 

To investigate the spectral statistics, the transport and 
the localization properties of the ERMH we apply the su- 
persymmetric technique (6|j^]. In the following, we choose 
a fixed realization of the diagonal elements rf R and re- 
strict the ensemble average to the random variables in £. 
This particular average is denoted by (• • -)f . Wc consider 
the generating functional 

F{ J) = H Dip exp [§<0(£ - H + (f + fe)A + J)tp] } . 

Here, ^ is a supervector with components ipj(R), which 
are themselves 8-dimensional supervectors with entries 



zi,z 1 ,xi,Xi, z 2,z 2 ,X2,X2, where the z v (x„) are com- 
plex bosonic (fermionic) variables. The diagonal matrix 
A has an equal number of eigenvalues +1 and —1 and de- 
scribes the grading into advanced and retarded Greens' 
functions. Furthermore, uj is a frequency, J a source ma- 
trix, and i/i is given by ip = ip'A. 

Choosing J appropriately fj],|| and taking derivatives 
of F(J) with respect to J at J = 0, one obtains ensem- 
ble averages of arbitrary products of Green's functions, 
in terms of which the above mentioned properties can be 
studied. Skipping most of the technical detail (for more 
information see we derive a nonlinear a model from 

the functional F(J). The physical implications are then 
discussed by comparing our formulation of the present 
problem with the supersymmetric description of a disor- 
dered metal given by Efetov . Performing the ensemble 
average, we get 

(exp(-i^))c = expf-i £ a(\R - R\) 



R,R 



istr[K (R)K(R)} 



(•5) 



where K(R) is a 8 x 8 super matrix given by K(R) = 
J2j v(j) ipj(R) $j{R)- The quartic term (||) can be de- 
coupled in the usual way by a Hubbard-Stratonovich 
transformation. This introduces a functional integra- 
tion over a field of 8 x 8 super matrices cr(R) with 
the same symmetries as K(R). Proceeding |J in anal- 
ogy to |Io| , [rT| , we apply a saddle point approximation, 
which is justified in the limit L\ 3> 1. Then we put 
a(R) — To + iTiQ(R), where Q(R) is an element of the 
orthogonal cr model space HQ and fulfills the nonlinear 
constraint Q(R) 2 = 1. The quantities Tq and Ti are 
determined by the implicit equation 



T + iT x A = -B v(j) J drj p (rj) 



x[E-r ] + ieA+-v{j)(T + iT 1 A)}- 1 , 
having the approximate solutions (in the limit L\ 1): 
T ~ -BqST 1 / dr) po(v)-Er- — > 

J E-7] 

Ti-wBoSpoiE) . 

Here, B = Er<\R\), S = E^O'), and V J dr, (_•_) 
denotes a principal value integral. From this and ^jfy 
we find the estimate T ± ~ U 2 /(W b L() < W b if U 
and W b are of the same order of magnitude. The ma- 
trix A R R = a(\R — R\) defines a d-dimensional gener- 
alization of a random band matrix Jl2[ with the typi- 
cal bandwidth L\. As a consequence of the Hubbard- 
Stratonovich transformation the coupling in the bilinear 
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term in the Q-field is given by the inverse, (A~ 1 ) RR . 
Therefore mainly the slow modes |^,[l2| (i.e. small mo- 
menta) contribute to the functional integral. Perform- 
ing a standard gradient expansion 

1@ and going over 
to the continuum limit, we can express the generating 
functional as F(J) = J DQ cxp[-£ 2 [Q] - AC(J)}. The 
effective action is given by 
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lV R Q(R)} 2 + f R (Q(R)) 
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+^0')( r o + *riQ)) , 



(0) 



where B 2 = l/(2d) J^r r2 a (\ R \) and A£(J) is the part 
of the action that accounts for the source matrix || . The 
Q-dependent "potential" Jr{Q) can be written in a more 
convenient form || 

f R {Q) ~ -i(7r/4) /i(r x /w) a>(£) str(QA) (7) 

with h(y) = iyY.\ 1 \<L c v i3)/[ l + njv(j)}, where L c « 
Li InLi is a cutoff length to be explained below. 

Before discussing this result, we mention that it is also 
straightforward to derive the ^-averaged local density of 
states from the functional F(J). It has the Breit-Wigner 
form 
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hr vV)]* + [hr lV (j)]* 



with an energy width Ti v(j) that depends on the relative 
coordinate j. The levels rf R of the product states \p1p2 > 
acquire a finite width (or inverse life time) Ti v(j) due 
to the interaction £. This width is of the order of Ti 
for \j\ < L\ and decreases exponentially like Y\ e~ 4 ^^ Ll 
for |j| 3> L\. This result demonstrates that product- 
states with 3> Li) have an exponentially large life 
time because both electrons are simply localized far away 
from each other without feeling the Hubbard interaction. 
Let us determine a critical length L c such that inside a 
volume of this diameter all pair states are well coupled. 
Equating the effective level spacing A e t f = Wb L~ 2d with 
the smallest possible level width Ti e ~ 2Lc / Ll , yields the 
estimate L c ~ L\\sxL\. Product states with \j\ > L c 
contribute to a discrete point spectrum while product 
states with \j\ < L c are well coupled and correspond to 
interaction-assisted pairs of the size L c . 

The dynamics of these well-coupled product-states is 
conveniently described in terms of the a model At 
this point we reiterate that our model describes diffusion 
(and localization) in the space of product states \p\p2)- 
The translation to coordinate space is straightforward 
provided we consider scales L > L\, where the \p1p2) 
are associated with well defined positions. The interpre- 
tation of our a model is greatly facilitated by the close 



similarities between (|^) and the a model for a disordered 
tal as derived by Efetov 0. Comparing (||) and (0) 



me 



with the standard a model in || , one can formally iden- 
tify an effective diffusion constant 2? e // and an effective 
local density of (diffusing) states v e f / , both of which de- 
pend on the frequency uj: 



I T\B 2 
tt B 2 



U 2 



v eff{u) = Po(E) h(T 1 /u>) 



(8) 
(9) 



We have introduced, via the Einstein relation, a formal 
"pair conductivity" a e f f , which does not depend on fre- 
quency. The function h(Tx/u>) can be interpreted as the 
number of states contributing to the diffusion. A detailed 
analysis Q yields the following limiting cases: 

*(ri/w)s , ri<H, 

/,(F, { [(L x /4) HT x /\u\)] d , r" 1 < H« Y x , 

[(L 1 /4)ln(r 1 T c )]' 1 , |u| < t- 1 . 

(10) 

Here, r c _1 is the effective level spacing A e // = 
Ti e - 2L c/ L i inside a blocks of size L c . The physical pic- 
ture of diffusing pairs makes sense in the regime \uj\ < Ti 
only. Furthermore, in formal analogy to the condition 
|cj| <C 1/r in a disordered metal Q (where r is the 
elastic scattering time), only those product-states with 
M Ti v(j) contribute to the diffusion. These remarks 
provide a physical interpretation of (^|) and (fio|). 

The analogy between the present problem and the 
problem of independent electrons in a disordered metal 
enables us to draw at least three important conclusions. 

First, the coupling constant f^e//2?e// = f' cr e// can 
be identified as a universal scaling parameter. The cor- 
responding scaling function is precisely the same as that 
of a disordered metal provided the latter is described 
by the "standard" a model ||. In particular, the per- 
turbative evaluation of the /3-function in 2 + e dimen- 
sions is equally valid for our present problem of 
diffusing or localized electron pairs (the term (Q) is not 
affected under the renormalization B). This first con- 
clusion also provides a rather rigorous justification for 
Imry's 
ture fll5 



application of the Thouless scaling block pic- 
For d = 1 we immediately recover Shepelyan- 
sky's original result jjj for the pair localization length 
L2, L% ~ o- e ff ~ (U /W 2 )L 2 . It is important to note 
that this result has been obtained by taking into account 
all, also the badly coupled (|j| ^> L\), pair states. In 
their study |Q Borgonovi and Shepelyansky argue that 
the badly coupled states should lead to a logarithmic cor- 
rection L 2 ~ Lf/lnLi. We cannot confirm this result. 

Second, we find that the pair dynamics is subdiffu- 
sive in agreement with recent numerical results of Bor- 
gonovi and Shepelyansky |^,^| who study two inter- 
acting kicked rotators. In the "pair-metallic" regime 
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(Li < L < L2 for d = 1,2) the a model can be 
treated perturbatively as in |J. The relevant diffu- 
sion propagator lZ(q,uj) = {<r e ffq 2 iuv e f /(a;)] -1 = 

[D e ff(uj) u e ff(uj) q 2 — iun/ e f f contains both the fre- 

quency dependent diffusion constant (||) and density 
of states (||). This uj dependence gives rise to sub- 
tle modifications of standard diffusion. Instead of try- 
ing to calculate the diffusion propagator lZ(R,t) = 
(27r)-( d+1 ) Jdiojdq Tl{q,uj)e^ qR -^ as a function of t 
and R we replace in a qualitative approximation u) by 
l/t. Therefore we expect the pairs to diffuse according 
to (see also @) 



(R 2 (t))~V eff (- t )t 



D ln(IV)- d i , Tf 1 «i <r c 
D ln(L!)- rf i , r c < t . 



We have used the notation D Q = (U 2 /W b ) h\- d . Obvi- 
ously, (R 2 (t)) increases weaker than linearly with time 
(subdiffusion) . The number of diffusing states given 
by the function h(Ti/uj) in (10) also depends on time. 
Putting h(Ti/ijj) = [L e f f (u})]^with L e ff(uj) the effective 
pair size, we get L e tf(l/t) ~ L\ ln(Fii) for T^ 1 < r c 
and L e ff(l/t) ~ L\ lnii for r c < i. This means that the 
pair size grows logarithmically with time as has also been 
found numerically in 0S. 

Third, calculating the two point correlation function 
Y"2(w), we can study the level correlations of the well- 
coupled product-states. The nearly localized pair states 
are mainly uncorrelated with the diffusing states and 
among themselves. Therefore their contribution to Y2(u>) 
essentially cancels out. For a finite system of size L the 

diffusive dynamics determines another energy scale, the 

(2) (2) 
"pair Thouless energy" E?>. In our case E c is given 

by the implicit equation V e ff(E^)/L 2 = E^. For 
lu < Ec the second term (0) of the action /^[Q] dom- 
inates the level correlations. The function Y2(oj) can be 
calculated in complete analogy with Q and we recover 
the random matrix result: 



r 2 M = y 2 (GOf! VAH) 



M « e, 



(2) 



Here, A(w) = [L d v e ff(uj)] 1 is the frequency dependent 
effective level spacing and Y2 GOE \r) is the universal 

spectral correlation function of the Gaussian orthogonal 

(2) 

ensemble. For higher frequencies, i.e. Ec <§C |w| <C Ti, 
also the first (kinetic) term of £2 [Q] has to be taken into 
account. The necessary perturbative evaluation of l2(w) 
proceeds analogously to the corresponding diagrammati- 
cal calculation of YiiijS) given by Altshulcr and Shklovskii 
HI. The result is 



Y 2 (u>) 



A 2 M 



co 2 \V eff (cu)/L 2 



d/2 



where V e ff(w)/L 2 can be interpreted as a frequency de- 
pendent Thouless energy setting the scale of the level 
correlation function. 

In conclusion, starting from an effective Hamiltonian, 
we have derived a nonlinear a model for two interacting 
electrons in a random potential in arbitrary dimension. 
Exploiting the analogy with Efctov's description of non- 
interacting electrons in a disordered metal, we identified 
a scaling parameter and investigated the level correla- 
tion function for the well coupled pairs. Furthermore, 
we analytically confirmed the numerical result that pair 
propagation is subdiffusive and that the pair size grows 
logarithmically with time. 
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